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^ , Abstract. We study the orbits of the various types of galaxies observed in the ESO Nearby Abell Cluster Survey. We combine 

the observed kinematics and projected distributions of galaxies of various types with an estimate of the mass density profile of 
[~^ ' the ensemble cluster to derive velocity-anisotropy profiles. Galaxies within and outside substructures are considered separately. 

Among the galaxies outside substructures we distinguish four classes, on the basis of their projected phase-space distributions. 
These classes are: the brightest ellipticals (with Mr < -22 + 5 log A), the other ellipticals together with the SO's, the early- 
^ I type spirals (Sa-Sb), and the late-type spirals and irregulars (Sbc-Irr) together with the emission-line galaxies (except those of 

early morphology). The mass profile was determined from the distribution and kinematics of the early-type (i.e. elliptical and 
SO) galaxies outside substructures; the latter were assumed to be on isotropic orbits, which is supported by the shape of their 
velocity distribution. The projected distribution and kinematics of the galaxies of other types are used to search for equilibrium 
^^ • solutions in the gravitational potential derived from the early-type galaxies. We apply the method described by Binney & 

^b ' Mamon as implemented by Solanes & Salvador-Sole to derive, to our knowledge for the first time, the velocity anisotropy 

f^ , profiles of all galaxy classes individually (except, of course, the early-type class). We check the validity of the solutions for 

"■^ ' P'ir) = [<Vj->(r)/<vJ^>(r)]''^, where <v^>(r)and<v^>(r)are the meansquared components of the radial and tangential velocity, 

"■^ ' respectively, by comparing the observed and predicted velocity-dispersion profiles. For the brightest ellipticals we are not able to 

construct equilibrium solutions. This is most likely the result of the formation history and the special location of these galaxies 
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Q ' at the centres of their clusters. For both the early and the late spirals, as well as for the galaxies in substructures, the data allow 

equilibrium solutions. The data for the early spirals are consistent with isotropic orbits (/8'(r) =1), although there is an apparent 

C/3 ■ radial anisotropy at ^ 0.45 ^200- For the late spirals an equilibrium solution with isotropic orbits is rejected by the data at the 
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> 99% confidence level. While fi'ir) » 1 within 0.7 ^200, P' increases linearly with radius to a value ^ 1 .8 at 1 .5 r2oo. Taken at 
face value, the data for the galaxies in substructures indicate that isotropic solutions are not acceptable, and tangential orbits 
are indicated. Even though the details of the tangential anisotropy remain to be determined, the general conclusion appears 
robust. We briefly discuss the possible implications of these velocity-anisotropy profiles for current ideas of the evolution and 
transformation of galaxies in clusters. 

Key words. Galaxies; clusters: general - Galaxies: kinematics and dynamics - Cosmology: observations 

1 . Introduction galaxies yield clues about the history of their accretion onto the 

cluster, and about the evolutionary relationships between them. 
The orbital characteristics of the various types of galaxies in 

present-day clusters can give unique information about the evo- The idea of shells of collapsing material around clusters 

lution of the clusters themselves, and about the formation and has been around since the work of Gunn & Gott ("1972"i. This 

evolution of their member galaxies. This is because clusters work stimulated several investigations, and indirect evidence 

are still accreting galaxies from their surroundings, and the de- for the infall of spirals into clusters has been accumulating over 

tails of this accretion process provide constraints for theories of the years. Moss & Dickens ri977 1 were the first to observe a 

cluster evolution. In addition, the orbits of the various types of diff'erence in the velocity dispersions of cluster early-type and 

late-type galaxies, followed by Sodre et al. JI989I I. and Biviano 

'^d offprint ^ests to: A. Biviano, e-mail: et al. tl992|. The difi^erent projected phase-space distributions 

bivianoOts . astro . it of early-type (red) and late-type (blue) galaxies was clearly es- 



* Based on observations collected at the European Southern tablished by Colless & Dunn J1996> and B iviano e t al. J1996> 
Observatory (La Silla, Chile) in the Coma cluster, while Carlberg et al. (■1997a") found it in 

** http://www.astrsp-mrs.fr/www/enacs.html the CNOC clusters at z ~ 0.3. The effect was studied in de- 
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tail for the clusters observed in the ESO Nearby Abell Cluster 
Survey (ENACS, hereafter; de Theije & Katgert [1999] paper 
VI; Biviano et al. 2002 paper XI), as well as in other clusters 
(see, e.g., Adami et al. .1998ai) . Galaxies with emission lines 
(ELG) provide a rather extreme example of the effect. The ELG 
are less centrally concentrated and have a higher dispersion of 
line-of-sight velocity than the galaxies without emission lines. 
This was first shown by Mohr et al. ( 1996 1 for the A576 cluster. 



and clearly demonstrated by Biviano et al. (119971 paper III) for 
the ENACS clusters. 

These results suggest mildly radial orbits of the late-type 
galaxies with emission lines, probably in combination with first 
approach to the central dense core. This interpretation would be 
consistent with the presence of the line-emitting gas which is 
unlikely to 'survive' when the galaxy crosses the cluster core. 
Indeed, Pryor & Geller fl984'l tried to constrain the orbits of 
Hl-deficient galaxies by noting that cluster-core crossing is a 



necessary condition for gas stripping, and Solanes et al. J2001> 
noted that the velocity-dispersion profile of Hl-deficient galax- 
ies is quite steep, suggestive of radial orbits. Support for the 
scenario of spiral infall into clusters comes from the the anal- 
yses of the Tully-Fisher distance-velocity diagram (TuUy & 
Shava[T984 Gavazzi et al. 1991 1. Indirect support comes from 
the numerical simulations that show that dark matter particles 
have a moderate radial velocity anisotropy, which increases out 
to the virial radius (e.g. Tormen et al. 1997 Ghigna et al. 1998] 
Diaferio ll999> . Radio or X-ray trails of cluster galaxies can 
also be used to constrain their orbits ('Merrifield ll998> . 

In the absence of full dynamical modelling, the analysis 
of the galaxy spatial distribution and kinematics can only sug- 
gest, but not really constrain, the nature of cluster galaxy orbits. 
This is because the projected spatial distribution, kinematics 
and mass model are coupled. So far, only a few full dynami- 
cal analyses of the orbital distribution of cluster galaxies exist. 
One reason for this is the relative paucity of detailed data on 
the kinematics and distributions of cluster galaxies, in partic- 
ular if several galaxy classes are considered. Another reason 
is that the orbital characteristics can only be inferred from the 
observed kinematics and distributions if the mass density pro- 
file of the cluster is known. The latter must be derived either 
from the distribution of light (with assumptions about the ra- 
dial variation of the mass-to-light ratio), or from the projected 
phase-space distribution of that subset of the galaxies for which 
the properties of their full phase-space distribution can be esti- 
mated independently. 

A first dynamical analysis of the orbits of cluster galax- 
ies was made for the Coma cluster by Kent & Gunn (1982"). 
Using several analytical mass models, these authors concluded 
that the galaxy orbits in the Coma cluster cannot be primarily 
radial, so that even at large radii a significant part of the ki- 
netic energy of the galaxies must be in the tangential direction. 
They noted that the range of the predicted velocity dispersions 
of the galaxies of different morphological types was only half 
that which is observed. Although a marginal result, this could 
indicate different distribution functions for the galaxies of dif- 
ferent types, and not just different energy distributions. Merritt 
J1987it used the same data to estimate the orbital anisotropy of 



the galaxies in the Coma cluster, for various assumptions about 
the radial dependence of the mass-to-light ratio. 

More recent dynamical modelling of galaxy clusters has 
led to the conclusion that the orbits of early-type galax- 
ies are quasi-isotropic, while those of late-type galaxies are 
moderately radial (e.g. Natarajan & Kneib 1996, Carlberg 
et al. 1997b. Mahdavi et al. il999 Biviano .2002. Lokas & 
Mamon l2003> . This picture is not supported by the analysis of 
Ramirez & de Souza 



I1998> who studied the deviations from 
Gaussianity of the overall distribution of the line-of-sight ve- 
locities of the galaxies. These authors concluded that the or- 
bits of ellipticals are close to radial, while spirals would have 
more isotropic orbits. However, van der Marel et al. (I 2000> 
and Biviano (p002 1 argue that the conclusion of Ramirez & 
de Souza is most likely due to erroneous assumptions in their 
modelling. 

One of the most extensive dynamical analyses so far was 
done for 14 'regular' galaxy clusters from the CNOC (Carlberg 
et al. ll997b1ll997c> . Adopting ad hoc functional forms for the 
3-D number density, the mean squared components of the ra- 
dial velocity, and the velocity anisotropy profile, Carlberg et al. 
(I1997bl I1997c> concluded that the velocity anisotropy is zero 
or at most mildly radial. The CNOC data were re-analysed by 
van der Marel et al. (2000 1, who used the method developed 
by van der Marel Cl994j . assuming a three-parameter family of 
mass-density profiles, and a set of constant values for the veloc- 
ity anisotropy, to determine the parameters in the mass-profile 
model from the best fit to the line-of-sight velocity dispersion 
profile. More recently, from the analysis of the projected phase- 
space distribution of ~ 15, 000 galaxies in the infall regions of 



eight nearby clusters (the CAIRNS project), Rines et al. (I2003> 
concluded that galaxy orbits are consistent with being isotropic 
within the virial radius. Note that neither van der Marel et al. 
(2000), nor Rines et al. C2003j distinguished among different 
cluster galaxy populations. 

In this paper we study the galaxy orbits in an ensemble 
cluster of 3056 galaxy members of 59 clusters observed in 
the ENACS. We use the 'inversion' of the Jeans equation of 
stellar dynamics, as derived by Binney & Mamon (1982), and 
we apply the solution method given by Solanes & Salvador- 
Sole ( 1990 hereafter S^). The analysis requires the mass profile 
M(< r), for which we use the estimate derived by Katgert et al. 
(2004 paper XII) for the same ensemble cluster Preliminary 
results were d iscusse d by Bivi ano et al. (tP^ EUU^ UJM\ . 
Mazure et al. (I2000> . Biviano (I2002> . and Biviano & Katgert 
(2003 ). 

In §|2lwe summarize the data that we use, describe the dif- 
ferent classes of cluster galaxies, and the construction of the 
'ensemble cluster' . In §|3we discuss the observational basis for 
our analysis, i.e. the projected and de-projected number den- 
sity profiles, and the velocity-dispersion profiles of the various 
galaxy classes. In §|4]we summarize the observed mass profile 
and the model fits that we used in the analysis of the orbits. In 
§|5lwe summarize the inversion procedure by which we derived 
the velocity-anisotropy profiles for the brightest ellipticals, the 
early spirals, the late spirals, and the galaxies in substructures. 
In §|6lwe describe the results of the analysis, and in §0we 
discuss the implications of the results for ideas about the evo- 
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lution of the clusters themselves, and about the formation and 
evolution of galaxies in clusters. Our summary and conclusions 
are given in § |S] Throughout this paper we use //q = lOO/i 
kms"' Mpc"'. 



2. The data, the galaxy classes, and the ensemble 
cluster 

Our analysis of the orbits of galaxies in rich clusters is based 
on data obtained in the context of the ENACS. Katgert et al. 
J 19961 []_998" papers I and V of this series, respectively) de- 
scribe the multi-object fiber spectroscopy with the 3.6-m tele- 
scope at La Silla, as well as the photometry of the 5634 galaxies 
in 107 rich, nearby (z $ 0.1) Abell clusters. After the spec- 
troscopic survey was done, a long-term programme of CCD- 
imaging with the Dutch 92- cm telescope at La Silla was car- 
ried out which has yielded photometrically calibrated images 
for 2295 ENACS galaxies. Thomas ( 2004 paper VIII) has used 
those images to derive morphological types, with which he also 
refined and recalibrated the galaxy classification based on the 
ENACS spectra, as carried out previously in paper VI. 

The morphological types derived by Thomas were supple- 
mented with morphological types from the literature, and those 
were combined with the recalibrated spectral types from paper 
VI into a single classification scheme. This has yielded galaxy 
types for 4884 ENACS galaxies, of which 56% are purely mor- 
phological, 35% are purely spectroscopic, and 6% are a combi- 
nation of both. The remaining 3% had an early morphological 
type (E or SO) but showed emission lines in the spectrum. With 
these galaxy types, Thomas & Katgert ( 2004 paper X) stud- 
ied the morphology-radius and morphology-density relations. 
These galaxy types also form the basis of the study of mor- 
phology and luminosity segregation (paper XI). 

In paper XI the galaxy classes were defined that must be 
distinguished because they have different phase-space distribu- 
tions. In particular, this applies to galaxies within and outside 
substructures. The membership of a given galaxy to a substruc- 
ture was determined using a slightly modified version of the 
test of Dressier & Shectman Cl988j . In this test, a quantity 6 
was computed for each galaxy, designed to indicate when the 
neighbourhood of the galaxy is characterized by a different av- 
erage velocity, and/or a smaller velocity dispersion than the 
cluster mean values (see paper XI for details). Galaxies with 
(5 < 1.8 were shown to have a very small probability of be- 
longing to substructures. On the other hand, only two thirds 
of the galaxies with (J > 1.8 really belong to substructures. In 
the present paper, we use 6 - 1.8 to separate galaxies within 
substructures from galaxies outside substructures. However, we 
also checked our results for the galaxies in substructures with 
6 > 2.2. Clearly, the 6 > 2.2 sample is smaller than the 5 > 1.8 
sample, but there is less contamination by galaxies outside sub- 
structures. The results for the 5 > 1.8 sample are confirmed 
from the 6 > 2.2 sample. Therefore, for the sake of simplic- 
ity, in the rest of this paper we only refer to 'galaxies in sub- 
structures' (or, more simply, 'Subs', in the following), meaning 
galaxies with <5 > 1.8, keeping in mind that the same results ap- 
ply for the galaxies with 6 > 2.2. 



In paper XI we showed that four classes of cluster galax- 
ies must be distinguished among the galaxies outside substruc- 
tures, on the basis of thei projected phase-space distributions. 
These are: (i) the brightest ellipticals (with Mr < -22+5 log/z), 
which we will refer to as 'E^j.', (ii) the other ellipticals together 
with the SO galaxies (to be referred to as 'Early'), (iii) the early 
spirals (Sa-Sb), which we will denote by 'Se', and (iv) the late 
spirals and irregulars (Sbc-Irr) together with the ELG (except 
those with early morphology), or 'Sj' for short. 

Summarizing, we consider 5 classes of cluster galaxies: 
E^j., Early, Se, Sj, and Subs, containing 34, 1129, 177, 328, 
and 686 galaxies, respectively. As explained in Appendix B.l 
of paper XII, corrections for incomplete azimuthal coverage 
in the spectroscopic observations and sampling incompleteness 
had to be applied in the construction of the number density pro- 
files. In order to keep these correction factors sufficiently small, 
galaxies located in poorly-sampled regions were not used and 
those have not been included in the numbers given above. 

The present analysis requires that data for several clusters 
are combined into an ensemble cluster, to yield sufficient sta- 
tistical weight. If clusters form a homologous set, the ensemble 
cluster effectively represents each of the clusters, provided that 
the correct scaling was applied. Support for the assumption of 
homology comes from the existence of a fundamental plane 
that relates some of the cluster global properties (Schaeffer et 
al. [T9931 Adami et al. 1998b paper IV; Lanzoni et al. l2004l . 
As shown by Beisbart et al. ( 2001 J . clusters with substructure 
deviate from that fundamental plane. Instead of eliminating all 
clusters with signs of substructure, we have chosen to consider 
separately those galaxies that are in substructures. 

As in papers XI and XII, we combined the data for 59 clus- 
ters with z < 0.1, each with at least 20 member galaxies with 
ENACS redshifts, and with galaxy types for at least 80% of 
the members (see Table A. 1 in paper XI). The resulting ensem- 
ble cluster contains 3056 member galaxies, for 2948 (or 96%) 
of which a galaxy type is available. The selection of cluster 
members was based on the method of den Hartog & Katgert 



(I1996> . and its application to the ensemble cluster is summa- 
rized in Appendix A of paper XII. We refer to papers XI and 
XII for details on the way in which the data for many clusters 
were combined. Those details concern the uniform method for 
the determination of cluster centres, and the correct scaling of 
projected distances from the cluster centres, R (with r2oo), and 
of relative line-of-sight velocities (with the global line-of-sight 
velocity dispersion o-p). The scaling with r2oo ensures that we 
avoid, as much as possible, mixing inner virialized cluster re- 
gions with external non-virialized cluster regions. Note that the 
scaling factors r2oo and o-p are computed using all cluster mem- 
bers. 

We assume that the ensemble cluster is spherically symmet- 
ric, not rotating, and in a steady state. As discussed at length in 
Appendix C of paper XII, these are reasonable assumptions for 
our ensemble cluster. 
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Fig. 1. Left: The best LOWESS estimate (solid line) of the projected number density l(K), within the 1-cr confidence interval 
determined from bootstrap resamplings (dotted lines), for each of the 5 galaxy classes, from top to bottom: E^j., Early, Se, Sj, 
Subs. The scale on the y-axis is arbitrary. Right: Same as left panel, but for the de-projected number density v(r). 



3. The number-density and velocity-dispersion 
profiles 

The observational basis for our study of the orbits of galaxies 
in clusters is provided by the projected number-density pro- 
files l(R), and the velocity-dispersion profiles crp{R) for the 5 
galaxy classes that we consider, viz. 'E|^j.', 'Early', 'Se', 'Sj', 
and 'Subs' (see §|2}. Here we summarize the steps involved in 
the determination of these profiles, and their de-projection. Full 
details can be found in Appendix B of paper XII. 

For the application of the Jeans equation - to derive the 
mass profile -, and its 'inversion' - to derive the velocity- 
anisotropy profile -, smooth estimates of number density pro- 
files, velocity-dispersion profiles and combinations thereof are 
required. We used the LOWESS technique (e.g. Gebhardt et al. 
I1994> to obtain smooth estimates of I(R) and o-p(R). Whereas 
Gebhardt et al. (1994 i appHed the LOWESS technique only to 
the estimation of a velocity dispersion profile, we also devel- 
oped a variant that produces a smooth estimate of the number 
density profile. 



The LOWESS technique yields estimates of I(R) and o-p(R) 
at the projected distance R of each galaxy. These estimates 
are based on a weighted linear fit to local estimates of pro- 
jected density and velocity dispersion. The linear fits typically 
involve between 30 and 80% of the data points, but with a 
weight that drops steeply away from the galaxy in question. 
The number density profiles, /(^)'s, were corrected for sam- 
pling incompleteness, assuming axial symmetry. Bootstrap re- 
samplings yield estimates of the 68% confidence limits (ap- 
proximately Icr-errors) of the LOWESS estimate. The pro- 
jected number density profiles I(R) of the 5 galaxy classes are 
shown in the left-hand panels of Fig. [2 together with their 68% 
confidence limits. 

In the Jeans equation as well as in its 'inversion' one also 
needs the de-projected 3-D number density v(r). In the right- 
hand panels of Fig. [Q we show the v(r)-profiles, as derived by 
de-projection via the Abel integral: 



v(r) 



1 r 



d/ AR 



AR V^TTTI 



(1) 
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This de-projection involves no assumptions other than spher- 
ical symmetry, the extrapolation of I(R) beyond the last mea- 
sured point towards large radii (for which we assume a tidal 
radius of 6.67 r2oo), and continuity of I(R) and its derivative 
at the last measured point. We checked that the de-projected 
profiles are essentially independent of the detailed form of the 
extrapolated I(R). 

In Fig.|2]we show the projected velocity dispersion profiles 
of the 5 galaxy classes as determined with the LOWESS tech- 
nique. In the same figure, we also show binned estimates of the 
velocity dispersion, where the value of crp{R) in each radial bin 
is computed using the robust biweight estimator (see Beers et 
al. [T99bl . 

To our knowledge, this is the first time that the number- 
density and velocity-dispersion profiles for these 5 cluster 
galaxy classes have been derived with such accuracy and in 
such detail. Therefore, we briefly comment on the qualitative 
nature of the different liR), v{r) and crp(R) before proceeding 
with the analysis. 

Among galaxies outside substructure, the E^j. have the 
steepest density profile in the centre, followed by the Early, the 
Se, and the Sj. This is a clear manifestation of the morphology- 
density relation (e.g. Dressler ll980> . and of luminosity segre- 
gation (e.g. Rood & Turnrose l 1 968l and paper XI). Interestingly, 
the density profiles of both Se and Sj decrease towards the 
cluster centre, a clear indication that these galaxies avoid the 
central cluster regions. On the contrary, E^j. are mostly found 
in the central cluster regions. 

The Subs galaxies have a number-density profile that is 
rather steep in the centre, but shows a weak 'plateau' at ~ 
0.6^200- Note that the number density profile of this galaxy 
class could, in principle, be biased by systematic effects due 
to the selection procedure of the members of substructures, 
which might result in a radius-dependent detection efficiency. 
A comparison of the de-projected number densities of the 
Subs-class galaxies and of the bulk of the galaxies outside sub- 
structures, viz. the Early-class galaxies (right-hand panels of 
Fig. 0, shows that, within ~ 0.6 r2oo, the two profiles have 
essentially identical logarithmic slopes. Beyond ~ 0.6 r2()o the 
number-density profile of the Subs galaxies is quite a bit flatter 
than that of the Early galaxies, until it steepens again beyond 
~ 1.0r2oo- This was already noted in paper XI. A comparison 
of the number-density profile of the Subs galaxies with that ob- 
tained by De Lucia et al. J2004t from their numerical models 
of substructures in cold dark matter haloes gives a similar re- 
sult. The logarithmic slope between 0.1 r2oo and 0.8 r2oo of the 
number-density of haloes with masses ~ lO^^M© is about -1.6, 
not very different from that of the Subs galaxies which is -1.5. 

The velocity dispersion of the E^j. strongly decreases to- 
wards the centre, with a slower but equally large decrease out- 
wards (remember that all velocity dispersions are normalized 
by the same, global velocity dispersion calculated for all galax- 
ies irrespective of type). The special formation history and lo- 
cation of the E|^j. at the bottom of the cluster potential well is 
reflected in their very low central velocity dispersion. In con- 
trast, galaxies of the Early class have a rather flat velocity- 
dispersion profile, changing by only » ±20% over the virial 
region. The velocity-dispersion profiles of Se andSj are rather 
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Fig. 2. The best LOWESS estimate (heavy fine) of o-p{R), to- 
gether with the 68% confidence levels (dashed lines), for each 
of the 5 galaxy classes, from top to bottom: E^j., Early, Se, Sj, 
Subs. The filled circles with error bars indicate binned biweight 
estimates of o-p{R). 



similar, starting at high values near the centre with a fairly rapid 
decrease out to r ^ 0.3 r2oo, and flattening towards larger pro- 
jected distances. Yet, the velocity dispersion of the Sj is larger 
than that of the Se (and, in fact, of any other class) at all radii. 
It is perhaps interesting to note that the velocity-dispersion pro- 
files of Se and Sj are remarkably similar to those of, respec- 
tively, the 'backsplash' and infalling populations of subhaloes 
found in the numerical simulations of Gill et al. (|2004i. 

Finally, the velocity-dispersion profile of the Subs class is 
very 'cold' and flat, even flatter and 'colder' than that of the 
Early class. One might wonder if this is due to the procedure 
by which the galaxies of the Subs class were selected, but it is 
very unlikely that the velocity dispersion of the Subs class is 
biased low by the selection. If anything, the actual velocity dis- 
persion of the subclusters is overestimated because the internal 
velocity dispersion of the subclusters has not been corrected 
for. In § 16. 41 we discuss several estimates for the real velocity- 
dispersion profile, i.e. corrected for internal velocity dispersion 
and possible bias due to the selection. 
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4. The mass profile 

In addition to the observed I{R)-, v(r)- and cr p(R)-pTofiles pre- 
sented in § 13 we also need an estimate of the mass profile 
M(< r) for a determination of the jS(r)-profiles. The mass pro- 
file that we will use here is the one that was derived in paper 
XII, from the number density and velocity-dispersion profiles 
of the Early-class galaxies. As discussed in detail in paper XII, 
the Early-class galaxies are likely to be in equilibrium with the 
cluster potential, as the formation of most of them probably 
antedates their entry into the cluster, so that they have had am- 
ple time to settle in the potential. In paper XII we also showed 
that galaxies of the Early class have a nearly isotropic veloc- 
ity distribution; this follows from an analysis of the shape of 
the distribution of their line-of-sight velocities. More specif- 
ically, assuming a constant velocity anisotropy for the Early- 
class galaxies, in paper XII we concluded that -0.6 $ /? $ 0.1, 
where 



Pir) = 1 



<V(>(r) 
<Vr>(r)' 



(2) 



and <Vj>(r), <v^>(r) are the mean squared components of the 
radial and tangential velocity (see, e.g., Binney & Tremaine 
11987 ). In this paper, we will often use the parameter y6' instead 
of yS to describe the velocity anisotropy, where /?' is defined as 
follows: 



li' = «y^^>l<v\>yi^ = (l-li)- 



(3) 



The constraint that we derived in paper XII for fi of the Early- 
class galaxies translates into/?' - l.O^gj^. 

For an isotropic velocity distribution (fi' - \ .0, or y6 = 0) 
the mass profile follows from the isotropic Jeans equation: 



r<v?> /dlnv dln<v?> 
M(<r)^ —\— — + 



dlnr 



where <Vj >(r) follows from: 



dlnr 



7rv(r) 



f 



2 1 p d[I(R) X fr2(R)] dR 

<v>(r) - ' 



dR 



V^^^72 



(4) 



(5) 



As with the de-projection of I{R), Eq.|5lrequires extrapola- 
tion of cTpiR) to the tidal radius (for details, see Appendix B.2 
in paper XII). 

The resulting M{< r), and its derivative p(r) are shown 
in Fig. 4 of paper XII. They are very well represented by 
a NFW profile (Navarro et al. I1997> with a scaling radius 

0.15 „ 
0.10 '^200- 



r, = 0.25+°'5 



5. The S^ method for the solution of /3'(r) 

Binney & Mamon < 19821 1 were the first to show that it is pos- 
sible to derive /?(r) when I{R), o-p(R) and M(< r) are known. 
S^ gave a practical recipe for application of the method, and we 
give a brief summary of their method to determine the velocity- 
anisotropy profile /3{r) for a given class of galaxies in equi- 
librium in a cluster gravitational potential with mass profile 
M(< r). Hereafter we give a brief summary of the S^ method 



(note that in this context we use /3 instead of/3' to be consistent 
with the earlier papers). 

The estimate of the mass profile M(< r) is used together 
with the estimate of the 3-D number density v(r) (derived 
from I{R) as before, see eq. [Q, to calculate ^(r) = -GM(< 
r)v(r)/r^. The observed functions o-p{R) and I{R) are used to 
derive H(R) - i I{R) crJ,(R), which in turn is used to calculate 
the function K{r) by the Abel integral: 



K{r) 



-r 



H(x) 



xdx 



VI^ 



(6) 



Using the functions ^(r) and K(r), one obtains the following 
two equations for <v^>(r) andy6(r): 



[3 - 2/3(r)] X <v2>(r) = ^ /^ **'(^) ^x - 

2 dK(r) 



(7) 



7Trv(r) dr 



and 



/3(r) <vXr) = 



-±^j;xMx)dx- 

, 1 dKjr) 



(8) 



;rrv{r) dr 
3K{r) , 3 



^J^K(x)dx 



7rnv(r) Krv 

from which <Vj^>(r) andy6(r) can be derived. 

The practical application of the method is far from trivial. 
First, one needs a smooth representation of the mass profile, 
which can be extrapolated confidently to large radii where we 
have not measured it. The extrapolation is done by using ana- 
lytic mass profiles that adequately fit the M(< r), such as the 
NFW profile (see paper XII). This ensures that the integral of 
^(r) in Eq.|S](whose upper integration limit we set to 6.67 r2oo; 
see §|3} is not problematic. Fortunately, ^(r) (which is nega- 
tive) asymptotically approaches with increasing r, and it does 
so with a sufficiently flat slope that the exact choice of the up- 
per integration limit and the analytic representation of M(< r) 
used for the extrapolation, do not influence the integral of *P(r) 
in a significant way. 

Secondly, one needs to extrapolate the observed velocity- 
dispersion profiles, without having very strong constraints. For 
each class, we check that dilTerent (plausible) extrapolations 
have no significant effect on the results of the S^ procedure 
within the observed radial range. 

A third important point is that Eq.|S]contains two integrals 
which have a lower integration limit of r = 0. Because it is 
quite difficult to determine the two integrands (r^*I'(r) and K{r)) 
at very small r from observations, a plausible interpolation of 
r^'^{r) and K{r) from the innermost measured 'point' to r = 
(for which both r^^(r) and K{r) are known from first princi- 
ples) is needed. We made a special effort to ensure plausible 
interpolations from the innermost point for which the data is 
available to r = 0, using low-order polynomials. 

It will not come as a surprise, given the equations involved, 
that it is practically impossible to give estimates of the for- 
mal errors in /3'(r) as derived with the S^ method. Approximate 
confidence levels on the /3'{r) of each galaxy class were there- 
fore determined by estimating the r.m.s. of four y6'(r), obtained 
by applying the S^ method to four subsamples, each half the 
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Fig. 3. The velocity anisotropy profile yS'(r) = 

<Vj>^^^/<Vj^>'^^, as derived for the galaxies of the Early 
class, using the mass profile that was derived assuming that 
the same galaxies havey6'(r) = 1. The shaded region indicates 
approximate 1-cr confidence levels, as obtained by considering 
subsamples half the size of the original sample. The value of 
jS'(r) is indeed quite close to 1, as it should. 

size of the original sample. The fact that each subsample only 
contains half the number of galaxies in the original sample, is 
likely to compensate for the fact that the four subsamples are 
not all mutually independent, which could lead to underestima- 
tion of the true confidence levels. 

We checked the robustness of our implementation of the S^ 
method as follows. We applied the S^ method to the galaxies of 
the Early class, adopting the mass profile that was determined 
using the same galaxies as isotropic tracers (see paper XII and 
§ 1^. Clearly, one should obtain fi{r) = 0, or, equivalently, 
/3'{r) = 1 (see eqs.|2and|3- The result is shown in Fig.|3 The 
shaded region indicates approximate 1-cr confidence levels, de- 
rived as described above. Indeed, we find a velocity anisotropy 
very close to zero with 0.85 < /3'(r) < 1.15 over the radial 
range < r/r2oo < 1-5. Deviation fromj6'(r) = 1 for the Early- 
class galaxies must be due to systematic errors arising from ex- 
trapolation uncertainties, and numerical noise in the inversion 
procedure (remember that our profiles are not analytic). Yet, 
the result in Fig.l^indicates that our implementation of the S^ 
'inversion' works quite well. 

We also applied a consistency test to all solutions that 
we obtained with the S^ method. I.e., we used the velocity- 
anisotropy profile /3{r) obtained by the S^ method for a given 
galaxy class, to determine the projected velocity dispersion 



profile through (see, e.g., van der Marel [T994l : 



v(r)<v>(r) 



and 



f 



v(^)M(<^) 



exp 



■r 



/3dx 

X 



R^\ vr<Vj^>(r)dr 



d^ (9) 



(10) 



I{R)(r,{R) - 2 1 - /3ir)- 

' Jr \ r^l Vr2 - R^ 

We then compared this predicted velocity-dispersion profile 
with the observed (Tp{R). In other words, we closed the loop. 



from observables and the mass profile to /3(r), then from /3{r) 
and the mass profile back to the observables. 

The observed and predicted o-p{R) are always in very good 
agreement (see §|6j, despite the fact that we cannot determine 
/3(r) beyond ~ 1 .5 r2oo, while knowledge of this function to very 
large radii is required to solve eq.^| The behaviour of /3(r) at 
large radii is not important since the number-density profiles 
of all galaxy classes drop sufficiently fast with radius. Even for 
the Sj, which have the shallower v(r), the effect of adopting two 
very different extrapolations of /3(r) to large radii (one derived 
from the analytical model proposed by Lokas & Mamon DOOll 
the other from the numerical simulations of Diaferio 19991 re- 
sults in a $ 10% variation at any point of the predicted crp(R). 

6. The velocity-anisotropy profiles 

We now investigate the orbits of the four classes of clus- 
ter galaxies that were not used to determine the mass pro- 
file, viz. E^j., Se, Sj, and Subs, in the gravitational poten- 
tial determined using the galaxies of the Early class. In other 
words: we try to construct equilibrium solutions for each of the 
galaxy classes, with physically acceptable velocity-anisotropy 
profiles. However, first we try to find solutions with isotropic 
orbits (or, /3' = 1). For this we need to solve eq. |9j using the 
v(r) of each class, and the mass profile M{< r) as determined 
using the Early-class galaxies, setting y8'(r) = 1. If the compari- 
son between the predicted and the observed velocity-dispersion 
profile yields an acceptable ;t'2, we conclude that the galaxies of 
the given class can be considered isotropic tracers of the cluster 
gravitational potential. 

After trying the isotropic solution, we then use the S^ 
method to solve for /3'(r). Note that, unlike Carlberg et al. 
( 1997^ ■1997c), van der Marel et al. ("200^), and Rines et al. 
(2003) we do not prescribe a functional form for j3'{r), nor do 
we assume a constant value for y6'(r). 

6.1. The brightest ellipticals 

The velocity-dispersion profile predicted for the E^j. class as- 
suming isotropic orbits is much flatter than the observed (Tp{R) 
(see Fig.|4j. We can reject the isotropic solution at > 99% con- 
fidence level (x^ = 98 on 4 data-points). 

Interestingly, abandoning the isotropy assumption does not 
help. I.e. there is no physical solution for which the E^j. are 
in equilibrium in the cluster gravitational potential (i.e. the S^ 
method predicts negative <v^> and /3{r) > 1 over most of 
the radial range covered by our observations). There are two 
straightforward interpretations of this result: either the galaxies 
of the E|jj. class are indeed out of dynamical equilibrium, or 
they do not fulfil the conditions for the application of the Jeans 
equation. We will return to this point in §0 

6.2. The early spirals 

For the galaxies of the Se class we do find acceptable equilib- 
rium solutions assuming an isotropic velocity distribution. The 
predicted velocity dispersion profile provides an acceptable fit 
to the observed crp(R) (x^ - 5.2 on 6 data-points, rejection 
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Fig. 4. The observed velocity-dispersion profile (Tp(R) of the 
E^j. galaxy class (dots with Icr errors), compared with the 
predicted o-p(R) (dashed line), assuming isotropic orbits in the 
gravitational potential determined from galaxies of the Early 
class. 



Fig. 6. The anisotropy profile, j0'(r) = [<v^>/<v^>]'^^ as de- 
rived for the galaxies of the Se class, through the S^ method. 
The shaded region indicates approximate 1-cr confidence lev- 
els, as obtained by considering subsamples half the size of the 
original sample. 
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Fig. 5. The observed velocity-dispersion profile o-p(R) of the 
Se galaxy class (dots with Icr errors), compared with the pre- 
dicted o-p{R) (dashed line), obtained by assuming isotropic or- 
bits in the gravitational potential determined from galaxies of 
the Early class, and with the predicted o-p(R) (solid line), ob- 
tained by using the velocity-anisotropy profile y6'(r) determined 
by the S^ method and shown in Fig.|5] 



probability of 61%). This profile is shown in Fig.|5] together 
with the observations and their Icr errors. Note that, although 
the data can be represented satisfactorily with isotropic orbits 
in the mass profile determined using the Early-class galaxies, 
the innermost values of crp{R) are somewhat underpredicted. 

The velocity-anisotropy profile of the Se class (determined 
via the S^ method) is shown as the solid line in Fig. |6l The 
velocity-anisotropy profile ;0'(r) is very close to unity near the 
centre, then rises to a maximum value of a: 1.8 at r a; 0.45 r2oo 
and then decreases again to reach a 1.1 at r/raoo ~ 1.5. As 
mentioned before, we checked the quality of this /3'{r) solution 



by calculating the implied velocity-dispersion profile, solving 
eq. |5]for this /3(r). The <Jp{R) predicted in this way from the 
yS'(r) indicated by the solid line in Fig.|6l is shown in Fig. 13 
also as a solid line. As expected, the latter is closer {x^ - 2.0 on 
6 data-points, rejection probability of 16%) to the observations 
than the isotropic solution (dashed line) but not significantly 
so, because the isotropic model akeady yields an acceptable 
fit to the data. As a matter of fact, the uncertainties on the y6'(r) 
profile determined via the S^ method are quite large, so that any 
deviation from the isotropic solution is not really significant. 

6.3. The late spirals+ELG 

For the galaxies of the Sj class we do not find acceptable 
equilibrium solutions assuming an isotropic velocity distribu- 
tion. This is illustrated in Fig. where the predicted crp{R) 
(dashed line) is clearly seen to provide a poor fit to the data 
(X^ = 18.2 on 6 data-points, rejection probability > 99%). 
Beyond R > 0.3 r2(K) the predicted velocity-dispersion profile 
is well below the observed values. Hence, purely isotropic or- 
bits are rejected. 

We then considered anisotropic solutions. The velocity- 
anisotropy profile of the Sj class (determined via the S^ 
method) is shown in Fig. |8] The profile is very close to unity 
out to r a 0.7r2oo, where it starts growing almost linearly with 
radius to reach a value of a 1.8 at r/raoo ~ 1-5. As usual, we 
checked the quality of the y6'(r)-solution by calculating the im- 
plied velocity-dispersion profile, solving eq.|5]for thisyS(r). The 
CTpiR) predicted in this way from the velocity-anisotropy pro- 
file indicated in Fig.|8] is shown in Fig. as a solid line. As 
expected, it reproduces quite well the observed o-p{R) of Sj. 

In the case of the Sj class the velocity-dispersion profile 
predicted with the yS'(r) obtained with the S^ method not only 
fits the data better than the isotropic case (this is also true for 
the Se-class galaxies), but it also does so in a significantly bet- 
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Fig. 7. The observed velocity-dispersion profile (Tp{R) of the 
Sj galaxy class (dots with \cr errors), compared with the pre- 
dicted o-p{R) (dashed line), obtained by assuming isotropic or- 
bits in the gravitational potential determined from galaxies of 
the Early class, and with the predicted crp(R) (solid line), ob- 
tained by using the velocity-anisotropy profile /3' (r) determined 
by the S^ method and shown in Fig.|S] 



ter manner (x^ - 2.6 on 6 data-points, rejection probability 
23%). Therefore, mild radial anisotropy is needed in order to 
put the Sj-class galaxies in dynamical equlibrium in the cluster 
potential. 

6.4. The galaxies in substructures 

As for the Sj -class galaxies, we do not find acceptable equilib- 
rium solutions for the galaxies of the Subs class if we assume 
an isotropic velocity distribution. As can be seen in Fig.|9l the 
predicted o-p{R) (dashed line) is way off the data (dots with 
error-bars; j^'^ - 118.5 on 6 data-points, rejection probability 
> 99%), and overestimates the observed velocity dispersion at 
essentially all radii. 

Using the S^ method for the Subs class,with the observed 
(Tp{R) we obtain the /?'(r) displayed as a solid line in Fig. 1101 
The orbits are tangentially anisotropic at all radii. As usual, 
we checked the fi'{r) solution in the space of observables; the 
predicted o-p(R) is in excellent agreement with the observed 
one (see Fig.|5J x^ - ^-^ on6 data-points, rejection probability 
60%). 

In the lower panel of Fig. ^| the shaded region indicates 
approximate 1-cr confidence interval, obtained as described be- 
fore. However, in this case the real confidence interval is prob- 
ably significantly larger, for two reasons. 

First, in using the observed velocity-dispersion profile of 
the galaxies in subclusters, we have ignored the internal veloc- 
ity dispersion of the subclusters. This means that the real veloc- 
ity dispersion is smaller than the observed one. We will make 
several assumptions for the (possibly R-dependent) value of the 
apparent internal velocity dispersion of the subclusters. In pa- 
per XI we estimated the internal velocity dispersion of the iden- 
tified subclusters, and obtained a value of ~ 400-500 km s \ 




Fig.8. The anisotropy profile, ;0'(r) = [<v^>/<v^>]'^^ as de- 
rived for the galaxies of the Sj class, through the S^ method. 
The shaded region indicates approximate 1-cr confidence lev- 
els, as obtained by considering subsamples half the size of the 
original sample. 




Fig. 9. The observed velocity-dispersion profile crp(R) of the 
Subs galaxy class (dots with Icr errors), compared with the pre- 
dicted a-p{R) (dashed line), obtained by assuming isotropic or- 
bits in the gravitational potential determined from galaxies of 
the Early class, and with the predicted o-p{R) (solid line), ob- 
tained by using the velocity-anisotropy profile yS'(r) determined 
by the S^ method and shown by the solid line in Fig.llOKlower 
panel). 



essentially independent of projected radius R. However, the 
true internal velocity dispersion of a subcluster is likely to be 
smaller, because the above estimate is biased high by galax- 
ies that do not belong to the subcluster but have been wrongly 
assigned to it by the selection algorithm. 

A more realistic estimate of a subcluster internal velocity 
dispersion is probably 250 km s"', a value close to the aver- 
age velocity dispersion of galaxy groups (see, e.g., Ramella et 
al. "1989"l. First, this constant value was subtracted in quadra- 
ture from the observed a-p{R) of the Subs class to produce the 
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Fig. 10. Upper panel: the observed velocity-dispersion profile 
o-p{R) of the galaxies in subclusters (solid line). The dashed 
line shows the result of deconvolving o-p{R) with an assumed 
constant internal velocity dispersion of 250 km s"'. The two 
other curves assume radius-dependent internal velocity disper- 
sions in the range 250^50 km s"' . Lower panel: the anisotropy 
profiles, fS'(r) = [<Vj>/<v^>]'^^ derived for the galaxies of 
the Subs class, through the S^ method with the observed I{R) 
and the four crp(7?)-estimates shown in the upper panel (with 
identical coding). The shaded region indicates approximate 1- 
cr confidence levels around the solution that uses the observed 
(Tp{R)-curve. These were obtained by considering subsamples 
half the size of the original sample. 



corrected velocity dispersion of subclusters shown as a dashed 
line in the upper panel of Fig.^| From this corrected crp(R), 
and using - as before - the observed liR), we derived the cor- 
rected version of y6'(r), indicated by the dashed Une in the lower 
panel of Fig. llOt . This second solution implies even stronger 
tangential anisotropy of the velocity distribution, which is not 
surprising since a larger fraction of the (smaller) line-of-sight 
velocities is required to balance the same cluster potential. 

However, it is possible that due to the selection procedure, 
or for physical reasons, we should not substract a constant 
value for the internal velocity dispersion of the subclusters. 
Therefore, we have assumed (rather arbitrarily) two different 



alternative solutions for the real velocity-dispersion profile of 
the Subs galaxies. These are shown as the dashed-dotted and 
dashed-triple-dotted curves in the upper panel of Fig.^| The 
former assumes a larger bias in the observed velocity disper- 
sion in the outer regions, while the latter mimics a larger bias 
in the central region. The important point in both assumptions 
is that the we must always deconconvolve the observed o-p{R) 
with at least 250 km s"' internal dispersion of the subclusters. 

For both assumptions about the real velocity-dispersion 
profile of the Subs galaxies, we calculated /3'{r), assuming - 
as before - that the observed I(R) is unbiased. The results are 
shown in the lower panel of Fig.^| where the same coding is 
used as in the upper panel. Not surprisingly, the evidence for 
tangential anisotropy of the Subs galaxies does not go away; if 
anything it gets stronger (in the most extreme cases, no physi- 
cal solution can be find beyond a certain radius). 

However, before we can accept this conclusion to be ro- 
bust, we must investigate the effect of possible biases in I{R). 
Without real modelling, we have considered two fairly extreme 
possibilities. In the first one we assumed that the plateau around 
0.6 r2oo in the observed v(r) of the Subs galaxies (see Fig.^and 
§ is (at least partly) an artefact due to the selection proce- 
dure. Consequently, we multiplied v(r) by a 'constant' factor 
of 2 below ~ 0.6 r2oo. The other possibility assumes that the 
logarithmic slope of v(r) below 0.6 r2oo is - 40 % flatter than 
actually observed. 

These two extreme assumptions about v(r) were combined 
with the various assumed estimates of (Tp{R) to solve for jS'(r) 
of the Subs galaxies. It appears that the conclusion of tangen- 
tial orbits is not affected by the different assumed shapes of v(r) 
below 0.6 r2oo, and is thus primarily driven by the low values 
of the velocity dispersion of the Subs galaxies. In other words, 
we find that, even without a detailed modelling of the selection 
effects and a possible dependence of the internal velocity dis- 
persion of the subclusters on radius, the conclusion about the 
tangential orbits is robust. The implications of this result will 
be discussed in §0 

7. Discussion 

Adopting the mass profile as determined from the galaxies of 
the Early class, we searched for equilibrium solutions for the 
other four classes. As a first step, we assumed isotropy of 
the velocity distribution, but subsequently we also solved for 
the anisotropy profile /3'(r) using the S^ set of equations (see 
§|5|l. For the E^j. class, we could not obtain equilibrium solu- 
tions, no matter what we assumed for their velocity-anisotropy 
profile. For the Se class the isotropic solution was found to 
be quite acceptable. Yet, the velocity-anisotropy profile of the 
Se class, as determined with the S^-method, shows a slight ra- 
dial anisotropy at r ^ 0.45 r2oo. For the Sj class, the isotropic 
solution is rejected. Their velocity-anisotropy profile, deter- 
mined with the S^ method, is close to zero out to r * 0.7 r2oo, 
and then increases almost linearly outwards, reaching a ra- 
dial anisotropy of/?' x 1.8 (corresponding to /3 ~ 0.7) at 
r X! 1.5 r2oo- For the Subs class the isotropic solution must also 
be rejected. Taken at face value, the data for this class imply 
substantial tangential anisotropy. However, this result may be 
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affected by systematic effects related to the selection of Subs 
galaxies. Until these effects have been modelled in detail, the 
conclusion of tangential anisotropy must be considered with 
some caution. 

Our conclusion that both early-type and late-type galaxies 
are in equilibrium in the cluster potential, with the latter on 
more radially-elongated orbits, is supported by several other 
studies in the literature. The larger velocity dispersion and/or 
the steeper velocity-dispersion profile of late-type galaxies with 
respect to early-type galaxies, have often been interpreted as 
evidence for infalling motions, and even for departure from 
virial equilibrium (Moss & Dickens 119771 Sodre et al ll989l 
paper III; Adami et al. ll998 al Solanes et al. 120011). However, 
Carlberg et al. J1997al l already pointed out that the latter does 
not need to be the case. They found that red and blue galaxies 
in the CNOC clusters are both in dynamical equilibrium in the 
cluster gravitational potential, but they were not able to con- 
strain the velocity anisotropy of these galaxies. 

From a more detailed analysis of the same dataset, van der 
Marel et al. (2000l were able to constrain the average veloc- 
ity anisotropy (assumed to be constant) of all CNOC cluster 
galaxies, to 0.75 < 13' < 1.2 (95.4% c.l.). This is a similar to 
what we found for the Early galaxies in paper XII, and from 
which we concluded that those have isotropic orbits. Mahdavi 



cluster centre can be understood with the model of Menci & 
Fusco-Femiano 



et al. J1999> showed that ELG in groups have an anisotropic 
velocity distribution, at the 95.4% c.l., with a best-fit constant 
P' ^ 1.8, whereas absorption-line galaxies in groups have a 
best-fit constant yS' ~ 1 .4, which is not, however, significantly 
different from unity. Both values seem somewhat higher than 
the values we find, which could indicate that the fraction of 
infalling galaxies is larger in groups than in clusters. 



Biviano (120021 . also using the ENACS dataset, concluded 
that, if absorption-line galaxies have zero anisotropy /?' = 1 .0, 
ELG have an average constant anisotropy of 1.3 < fi' < 1.6 
(68% C.L). This range is in reasonable agreement with our 
result for the velocity anisotropy profile of the Sj galaxies 
(see Fig. |8j, considering that many Sj galaxies are found at 
large radii, where their radial velocity-anisotropy is largest. 
Natarajan & Kneib (1996 1 concluded that galaxies in A2218 
have tangential orbits in the central region, and radial outside, 
with an anisotropy profile resembling the one we find for Sj. 
Finally, Ramirez & de Souza (1998'l and Ramirez et al. (2000) 
concluded that early-type galaxies have more eccentric orbits 
than late-type galaxies, but their result arises from incorrect as- 
sumptions in their method, as discussed by van der Marel et al. 
( 12000) and Biviano (2002). 

The present analysis is the first to consider the orbits of 
5 distinct classes of cluster galaxies. Using the Early class as 
a reference, we find that 3 of the remaining 4 are in dynam- 
ical equilibrium within the cluster gravitational potential; this 
is manifestly not the case for the E^j.. The most likely expla- 
nation for our failure to find solutions of the collisionless Jeans 
equation for the Ei^j., is that the Ej^j. either formed very near the 
cluster centre, or moved there by losing kinetic energy subject 
to dynamical friction. At the same time, they probably have 
grown through merging with other galaxies. These processes 
lead to a loss of the orbital energy of these galaxies. As a mat- 
ter of fact, the very low velocity dispersion of the E^jj. at the 



||1996> . which is a solution of the collisional 
Boltzmann-Liouville equation, and hence accounts for galaxy 
collisions and merging processes. 

The conclusion that the Early-class galaxies have a nearly 
isotropic velocity distribution is not surprising, given the large 
body of evidence indicating that ellipticals are an old cluster 
component. If they form and become part of the cluster be- 
fore it virializes, they can obtain isotropic orbits through vi- 
olent relaxation. From the distribution of the ratio rpgri/fapo 
of the dark matter halos in their simulations of rich clusters. 



Ghigna et al. ( 1998 ) concluded that about 25% of the halos are 
on orbits more radial than 1:10, where the median ratio is 1:6. 
Comparison with our result is not immediate, but of the galax- 
ies outside substructures 36% belong to the Se and Sj classes. 
Since about two-thirds of those are late spirals or ELG, which 
are the galaxies showing most of the velocity anisotropy, these 
could indeed correspond to the halos with orbits more radial 
than 1:10 in the rp(.„7rapo -ratio. 

The increase of the radial-velocity anisotropy with radius of 
the Sj (see Fig. |8} is a feature commonly found in numerical 
simulations of dark matter haloes. E.g., the numerical simula- 
tions of Tormen et al. ( I1997> predict an increasing radial veloc- 
ity anisotropy from r/r2()o ~ 0.3 outwards, reaching yS' ~ 1.8 
at r/r2oo ~ 1-5, and the numerical simulations described by 
Diaferio et al. |2001 ) predict a similar, though somewhat more 
irregular, behaviour of the velocity anisotropy profile, with a 
maximum anisotropy of /3' ~ 1 .4. This anisotropy profile re- 
sults from infall motions of the field haloes into the cluster, and 
from the subsequent isotropization of the velocity distribution 
of these haloes as they move towards the denser cluster centre. 

The similarity of the /3'(r) of the Sj-class galaxies and of 
the dark matter haloes in the models is quite interesting and it 
probably means that the Sj galaxies still retain memory of their 
infall motion from the field. The fact that a large fraction of the 
Sj have emission-lines indicates that they have not yet lost their 
gas as a consequence of tidal stripping, galaxy collisions, or 
ram pressure. Hence it is unlikely that the Sj we observe have 
spent much time in the hostile cluster environment, and many 
of them could indeed even be on their first cluster crossing. 

On the other hand, the small velocity anisotropy of Sj near 
the centre probably reflects the fact that the galaxies we iden- 
tify as Sj must avoid, or have avoided, the central region. Those 
Sj that have a significant radial anisotropy near the centre will 
cross the very dense central cluster regions, where they can- 
not survive and get disrupted, either to form dwarfs, or to con- 
tribute to the diffuse infra-cluster light (Moore et al. ll999al . As 
a matter of fact, Sj are not found in the cluster central regions 
(see Fig. [2 and paper XI). The existence of faint spiral struc- 
tures in some dwarf spheroidals has now been demonstrated 
(Jerjen et al. l2000l Barazza et al.'2002' Graham et al.'200T). 
Interestingly, Conselice et al. (2001 ) found that, in the Virgo 
cluster, the dwarf spheroidals have a velocity distribution more 
similar to that of the spirals, than to that of the ellipticals. 

Our analysis shows that an acceptable equilibrium solution 
exists for Se with zero velocity anisotropy. Hence, it is likely 
that these galaxies are not very recent arrivals, since there is 
no evidence for memory of their initial infall motion. These 
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Se galaxies are more likely to survive the hostile cluster envi- 
ronment than Sj, because of their higher surface brightness (see 
paper X; Moore et al. 1999a i. This is consistent with the results 
of numerical simulations showing that clusters contain red disk 
galaxies that, after accretion from the field, attain dynamical 
equilibrium in ~ 1-2 Gyr (Diaferio et al. 2001 1. 

Additional indirect support for the scenario described 
above comes from the similarity of the Se and the Sj velocity 
dispersion profiles with those of, respectively, the 'backsplash' 
subhaloes, and the subhaloes on first infall, identified by Gill et 
al. (2004 1 in their N-body simulations of galaxy clusters. This 
similarity suggests that many Sj could be on first infall, while 
the Se at large radii have already crossed the cluster core. 

In paper X it was argued that Se are likely to be the progen- 
itors of SOs. This conclusion is based on three different pieces 
of evidence: (1) the strong increase of SOs in clusters since 
z ~ 0.5 (Dressier et al. 1997; Fasano et al. 2000 1, accompanied 
by a similar decrease of the spiral fraction; (2) the morphology- 



density relation (Thomas et al. 120041 hereafter paper IX) which 
shows that the local projected density around Se is smaller 
than around SOs; and, (3) the strong similarity of the bulge lu- 
minosity of Se and SOs (paper X). If Se transform into SOs 
and if the velocity distribution of SOs is isotropic (see paper 
XII), it is only natural that the velocity distribution of Se is 
also isotropic. Otherwise, the timescale of the morphological- 
transformation process should be similar to that of the velocity 
isotropization. 

Even if the isotropic solution is perfectly acceptable for the 
Se, the data, when taken at face value, imply some radial ve- 
locity anisotropy. Although the significance of the anisotropy is 
rather low, we are tempted to speculate about a possible cause, 
if the anisotropy at r/r20() ~ 0.45 were real. Galaxies with ra- 
dial velocity-anisotropy, moving on radially elongated orbits, 
will move relatively fast near the cluster centre. It is possible, if 
not likely, that the high radial velocity is a necessary condition 
for Se to avoid impulsive encounters and thus transformation 
into SOs in the central high-density cluster region. 

Finally, galaxies in substructures provide an intrigueing 
view into the processes that are important in the formation of 
clusters. Recently, several groups have studied the properties 
of substructures within dark-matter haloes over a range of to- 
tal masses that includes those of rich clusters (e.g. De Lucia 
et al. 2004; Taylor & Babu F2004j . We already mentioned the 
agreement between our radial number-density profile of Subs 
galaxies and the radial distribution of 10^^ Mq substructures in 
the models of De Lucia et al. ( 2004 see § |3}- Those models 
also show that the more massive substructures are preferen- 
tially located in the external regions of their parent haloes. This 
is most likely due to tidal truncation and stripping of substruc- 
tures that reach the dense central regions. In addition, orbital 
decay can also contribute to this mass segregation (e.g. Tormen 
et al. [1998.1. The apparent paucity of Subs galaxies in the in- 
ner regions of our clusters may thus well be the result of mass 
segregation, instead of selection bias. 

Taylor & Babul (2004i discuss the evolution of the orbits 
of the infalling substructures, and they conclude that disrup- 
tion occurs sooner for more radial orbits. This will lead to a 
tangentially anisotropic distribution of orbits of the surviving 



substructures, which is exactly what we find. So, even if the de- 
tails of the tangential anisotropy of the Subs galaxies requires 
additional modelling, the result itself appears robust and not 
unexpected or implausible. 

8. Summary and conclusions 

We determined the equilibrium solutions for galaxies of the 
4 classes that were not used as tracers of the cluster poten- 
tial. For this, we solved the inverse Jeans equation, using the 
method of S^. We found equilibrium solutions for galaxies of 
the Se, the Sj, as well as the Subs classes, but not for galaxies 
of the E|^j.-class. The equilibrium solution found for galaxies 
of the Se class was found to be consistent with them being on 
isotropic orbits, except perhaps just outside the cluster central 
region. On the other hand, isotropic solutions were found not to 
be acceptable for galaxies of either the Sj or the Subs classes. 
Galaxies of the Sj class were found to be on mildly radial or- 
bits, with the radial velocity-anisotropy increasing outwards. 
On the contrary, tangential orbits seem to characterize galaxies 
of the Subs class, but the significance of this result is difficult 
to assess in view of possible systematics effects we have con- 
sidered. 

Our results support hierarchical models for the build-up of 
galaxy clusters (see also paper XII). Our results also constrain 
the evolutionary history of cluster galaxies. They are consistent 
with, if not suggestive of, a scenario where the very bright ellip- 
ticals form very early, and sink to the bottom of the still form- 
ing cluster potential well, losing orbital energy. In our scenario, 
the less bright ellipticals, together with the SOs (the Early-class 
galaxies), were already part of the cluster at the epoch of its 
formation, and developed isotropic orbits through the process 
of violent relaxation, or have lived sufficiently long in the clus- 
ter to have lost any memory of original radial infall motions, 
through isotropization of their orbits. 

This is probably also the case for the early spirals, which 
make them acceptable candidates for being the progenitors of 
SOs, also in view of their structural properties. We speculate 
that some early spirals near the cluster centre have managed 
to escape transformation into SOs as a result of a selection ef- 
fect in the velocity distribution. Finally, many late spirals and 
emission-line galaxies (excluding those of early morphology) 
are likely to be field galaxies recently arrived into the cluster. 
Their radial infall motions are gradually isotropized as they 
approach the cluster centre, until they get disrupted or trans- 
formed into dwarf spheroidals as a consequence of collisions 
and, in particular, tidal effects. 

The galaxies in substructures apparently avoid the central 
regions and they appear to be on tangential orbits. Although 
some modelling remains to be done to assess the details of 
the implied anisotropy profile, the conclusion of tangential 
anisotropy appears to be robust. Interestingly, both effects are 
also seen in numerical simulations, and they result from the 
mechanisms that 'destroy' the substructures as they get nearer 
to the cluster cores. 
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